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Abstract. 

The Dirac mass-vortex at the surface of a topological insulator or in graphene 
is considered. Within the linear approximation for the vortex amplitude's radial 
dependence, the spectrum is a series of degenerate bound states, which can be classified 
by a set of accidental SU{2) and supersymmetry generators (Herbut I F and Lu 
C K 2011 Phys. Rev. B 83 125412). Here we discuss further the properties and 
manifestations of the supersymmetry of the vortex Hamiltonian, and point out some 
interesting analogies with the Runge-Lenz-Pauli vector in the non-relativistic hydrogen 
atom. Symmetry breaking effects due to a finite chemical potential, and the Zeeman 
field are also analyzed. We find that a residual accidental degeneracy remains only 
in the special case of equal magnitudes of both terms, whereas otherwise it becomes 
removed entirely. 
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1. Introduction 

Accidental degeneracy has been an interesting issue in mathematical physics, with 
the classic examples ranging from non-relativistic hydrogen atom to multi-dimensional 
harmonic oscillator. The fact that states labeled by different quantum numbers share 
the same energy is explained by identifying hidden constants of motion, unrelated to 
spatial symmetries that may be present in the problem. In hydrogen atom, for example, 
these are the three Runge-Lenz-Pauli operators which, together with the usual angular 
momentum, constitute a closed so(3) x so(3) Lie algebra, and so imply the degeneracy 
among all the levels of a given principle quantum number. [H |2] For the hydrogen spectra 
in the relativistic limit, however, the linear dispersion in the Dirac Hamiltonian together 
with its matrix structure complicate the degeneracy pattern and hinder the search for 
the hidden accidental constants of motion. [31 H] 

Recently, the interface between topologically trivial and non-trivial insulator has 
become a topic of many investigations, because of its hosting of the quasi-relativistic 
chiral Dirac fermions.[5] Depositing a type-II superconductor onto the surface of a 
topological insulator, for example, could by the proximity effect induce a mass-term 
for the Dirac fermions, proportional to the superconducting order parameter. [6] It 
is well known that, in general, in the vortex configuration of such a mass-term, 
there exists a special topologically protected state in the Dirac spectrum at zero 
energy. [7] Many authors have recently elaborated on the properties of such zero-modes 
in condensed matter systems, and derived several unexpected consequences of their 
existence. [El El [TOl HI [El [131 El El [16] 

Obtaining the entire bound spectrum for a general vortex profile A(r), however, 
is typically not possible analytically, and one has to resort to numerical methods. [17] 
Nevertheless, the corresponding Dirac Hamiltonian becomes analytically solvable once 
the amplitude of the vortex is linearized near the origin. [1H1I19] In our previous paper. [19] 
we showed that due to the presence of a closed Clifford algebra associated with Dirac 
matrices the energy spectrum even in this "harmonic approximation" is non-trivial, 
and, in particular, exhibits degeneracies due to the hidden and accidental SU{2) and 
supersymmetries. The eigenvalues are given hy E = ±u\/N, with an energy scale 
u, and with integer corresponding to the total number of fictitious fermions and 
bosons that can be formally defined in the problem, and then distributed over two 
(in two dimensions) quantum states in all possible ways. The associated degeneracy 
then comes out to be 2A^, and when > 1 there are two distinct multiplets of the 
effective "angular momentum", with j> = A^/2 and j< = A^/2 — 1, reminiscent of the 
nonrelativistic spectrum of the hydrogen atom, i2j Such a bound-state spectrum could 
be observable by the scanning tunneling microscopy, and would determine the specific 
heat of the vortex, for example. 

In this paper, we investigate further the properties and the algebra of the 
supersymmetry generators, which are responsible for the degeneracy between different 
multiplets of the angular momentum. These generators commute with the square of the 
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Dirac vortex Hamiltonian, and their emergence in the problem is closely related to the 
existence of different square roots of the identity matrix. We also establish a partial 
analogy with the Runge-Lenz-Pauli vector for the text-book hydrogen atom. [Tl|2] Besides 
its relevance to topological insulators, our study is also pertinent to the insulating or 
superconducting graphene, where the low-energy excitations are also Dirac fermions, and 
the mass-gap can in principle be generated dynamically, as well as externally . [TOl [TTl [H] 

In addition to the usual symmetry of the Dirac spectrum around zero energy, the 
chiral symmetry of the Dirac Hamiltonian is also responsible for the appearance of the 
linear combination of odd and even occupation numbers of the fictitious fermions in the 
bound-state wave functions. We use this property to investigate the symmetry breaking 
effects of the chemical potential and the Zeeman magnetic field, which in reality are 
present both in graphene and in topological insulators. We find that the degeneracy is 
completely lifted in general, but some residual degeneracy remains in the special case 
when the two contributions are equal in strength. 

The paper has the following organization. In Sec. Ej the supersymmetry generators 
and their superalgebra are introduced in terms of the associated supercharge. Sec. [3] is 
devoted to explicitly constructing these symmetry generators. A comparison with the 
Runge-Lenz-Pauli vector is also provided. The effects of the finite chemical potential 
and the Zeeman field, and their connection with the chiral operator are investigated 
in Sec. HI Finally, the discussion of validity of the linear approximation made for the 
vortex amplitude, and comments on related problems in the literature are offered in Sec. 
[5l Some explicit wave functions are exhibited in Sec. El 

2. Emergent supersymmetry from H'^ 

We define the vortex configuration 



in the two-component order parameter A(r), which when uniform provides a mass-gap 
for the two-dimensional Dirac fermions. Such an order parameter may in principle 
arise from the proximity-induced superconductivity in both graphene and topological 
insulator, for example, (r, 6) are the usual polar coordinates in the plane. The 
winding number of the vortex is assumed to be unity. The spatial profile of the the 
order parameter's amplitude will be taken to be qualitatively similar to the simple 
|A(r)| = Aqo tanh (r/^), which features the generic linear dependence on the radial 
coordinate near the origin and the saturation to a finite value A^o at infinity. ^ is the 
characteristic (coherence) length. 

With the above qualifications in mind, the linearized Dirac Hamiltonian with the 
mass-vortex may be written as 




(1) 




(2) 
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where the set of four 4x4 Hermitian matrix i = 1,2, constitutes the 4- 

dimensional Chfford algebra. The fifth matrix F = aia2/3i/32 is responsible for the 
chiral symmetry, since it evidently anticommutes with the Hamiltonian: {F, H} = 0. 

A supersymmetric representation of the Dirac Hamiltonian in Eq. (|2]) is possible 
both because of the existence of the standard bosonic representation for the operators 
of coordinates and momenta. 




&. = y^a:. + zy^p., (3) 



and of the fermionic representation for the four Dirac matrices, 



(4) 



2 

Here = {ai,a]} = 6ij, [bi,bj] = {ai,aj} = 0, as usual. The two transformations 

together yield the following supersymmetric Hamiltonian, 

H = uj{a\bi + alb2 + aib\ + a2bl) , (5) 



where the energy u = y2cAoo/^- Hereafter, oo will be assumed to be unity, and will 
not be shown exphcitly unless needed. The energy eigenvalues in this representation 
become obvious, because 

H'= J2{ala, + blbi), (6) 

i=l,2 

turns out to be a simple sum of fermion and boson particle numbers. Therefore, the 
energy levels can be labeled by the total particle number A^, and various ways of 
distributing these fictitious particles according to their statistics leads to the degeneracy 
of 2iV.[19] It is also easy to see that the zero-mode in this representation corresponds 
to the nondegenerate vacuum |0), which is also an eigenstate of the chiral symmetry 
operator F with the eigenvalue of unity. 

In order to see how the underlying supersymmetry emerges, we first note the 
commutators [H^, aj] = aj and [H^, b]] = 6| for z = 1,2. We then search for 
other operators whose square has the identical commutation relations with the particle 
operators. With the help of identity [A^^X] = {A,[A,X]} = [A,{A,X}] for either 
X = b\ or X = aj, one readily finds the following three such additional supersymmetry 
operators: 

Ai= aijk{a]bk + b]ak), (7) 

j,k=l,2 

which satisfy A"^ = (without the summation convention). In the above, the index i 
denotes the i-th Pauli matrix, whereas the subscripts (j, k) label the matrix elements. 
The existence of these operators is linked to the fact that the two-by-two identity 
matrix I2 has four distinct square roots. Together with H, the four supersymmetry 
operators may be decomposed into the supercharges and their Hermitian conjugates. pO] 
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{H,A) = Qfi + Qj^- The time-component /x = corresponds to H and the rest of 
components refers to A. The supercharge is defined as 

Q^l= J2 b]i'^f^)jkO'k- (8) 

We now determine the algebra satisfied by the supersymmetry operators. First, due 
to the quantum statistics associated with the a's and 6's, the supercharges are mutually 
ant icommut ing : 

{Q^.,Q.} = 0. (9) 

Next, it is useful to invoke the anticommutator 

{albj, bjam} = ajam^ji + h\hj5im , (10) 

between boson-fermion mixed operators to compute the anticommutators between Q's 
and 

{Qti^Ql] = {h]{a^ay)jkhk + a]{ai,a^)jkak) , (11) 

which will be crucial in the subsequent generation of the constants of motion. Eq. (fTTj) 
leads to the following simple superalgebra for the components of the operator A: 

{Ai,A,} = 2H^5^j. (12) 

When the time- component H is involved in Eq. (fTT]) . on the other hand, one obtains 
the generators of the SU (2) symmetry of the Dirac Hamiltonian as 

{H,A} = U, (13) 

defined as.[T9] 

'^i = ^ (^i,jk{a]ak + b]bk) . (14) 

^ j,k=l,2 

From Eq. f|T2|) and f|T3|) . we then see that J is a constant of motion, since 

[J„H] = ^[{H,A,},H] = 0. (15) 

Next, the algebra between J's and A's may be displayed in terms of commutators, 
which also make the vector nature of the latter under the SU{2) symmetry manifest. 
The following identity, 

[XjXj, XjXm] = XjXmSjl — XiXjSim , (16) 

for both fermion X = a and boson X = b can be used to derive the standard SU{2) 
commutator algebra between the angular momentum generators, 

[Jj, Jj] = ieijkJk ■ (17) 

From the commutator between distinct species, 

[XlY„Y^Y^]=XlY^5,i, (18) 
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for X = a{h) and Y = b{a), on the other hand, one can show that A is indeed a vector 
under the SU{2) algebra, i. e. 

[Ai, Jj] = ieijkAk . (19) 

Eqs. (fT2!) . (fT3|) . (fT7|) and (fT9l) constitute the superalgebra that is behind the 
degeneracy of the spectrum of H. The reader should note that Eq. (fT3|) . (fT7|) and 
(fT9l) are mutually consistent: this can be seen, for example, through replacing with 
Ji = {H,Ai}/4 to show that [Ji, J2] = {[Ai, J2],i/}/4 = iJg. In order to include the 
chiral symmetry operator T and complete the algebra, we note that {F, at} = [T, bi] = 0, 
which yields the final relations. 



{T,H} = {T,A} = 0, (20) 

[r, J] = . (21) 

Let us remark in passing that the above derivation can be extended to higher 
dimensions, to consider, for example, the three-component mass hedgehog in the three- 
dimensional Dirac equation. [T6| [T9| [2T| This case would involve distributing fermions 
and bosons over three quantum states to obtain the spectral degeneracies. The essential 
modification would then be the replacement of the matrices {I2, <?} which appear in 
Eq. ([7]) and (fT4|) with another set, consisting of an identity matrix and the eight Gell- 
Mann matrices. [22] Again, the equalities Eq. (ITOl) . (|T6l) . and ( ITSl) are still valid between 
three boson and three fermion operators. Therefore, the commutation relations in Eq. 
( fTTl) and ( fT9l) still hold, except that the structure factor associated with SU{3) algebra 
should be involved. Although [Af,!!"^] = also holds here, Af does not coincide with 
as in SU{2) case. Nevertheless, the square of Hamiltonian still can be expressed as 
the sum = (3/32) ELi 4- 



3. Representation of supersymmetry operators and connection with 
Runge-Lenz-Pauli vector 

The degeneracy pattern of the spectrum of the Hamiltonian in Eq. (j2]) (Fig.l in Ref.|19]) 
is similar to that of the bound spectrum in hydrogen. It is therefore of interest to 
investigate how the states within each eigenspace with E = \fN are connected to each 
other via the application of the operators J and A. In addition, since the subspace of 
the second excited level with N = 2 has a one-to-one correspondence with that of the 
first excited level with = 2 in the hydrogen atom, we should be able to see the exact 
correspondence in algebraic relations between the two cases. 

First, a matrix representation for vector operator A of Eq. ([7]) will be constructed. 
We work in the basis of eigenstates \E,j,m) labelled by the respective eigenvalues of 
H, and J3. We recall that for > 1 there are two allowed magnitudes of angular 
momentum j> = A^/2 and j< = j> — 1 within the degenerate subspace with E = y/N . [T^] 
The resultant degeneracy is then 2A^. The chiral operator F also connects the states 
with opposite signs of the energy, as usual. 
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According to Eq. (1191) . commutes with J3 and H^. When > 1, it suffices then 
to consider the subset consisting of the four states {\±E,j^,m), \±E, j^,m)}. We now 
show that: 



A3|E,j>,m) = — |E,j>,m) + - -^1 - ^,J<,m) . (22) 
First, bracketing {ff, ^3} = 4 J3 between a pair of states leads to 

2777, 

(E,j>,m|A3|^,j>,m) = — , (23) 

and to the conclusion that A^\E, j^^m) is orthogonal to |i?,j<,m). Next, due to Eq. 
( 120|) . the matrix element 

(-E,j>,m|A3|E,j>,m) = -(-E,j>,m|A3|E,j>,m)* (24) 

must be purely imaginary, or zero. Since ^3 in Eq. ([7]) has only real coefficients, one 
concludes that y43|i?,j>,m) is orthogonal to | — E,jy,m) as well. Consequently, 

A3IE, j>,m) = —\E,j>,m) + d\-E,j^,m) (25) 
with a yet undetermined coefficient d. Similarly, 

2777- 

A3I -E,J<, 777) = - — I -E,J<, 777) + C?|E,J>,777). (26) 

Using then A"^ = fixes the value of the coefficient d to d"^ = E"^ — Am'^/E'^, as in Eq. 
(1221) . Utilizing the commutators [^3, J±\ = ±y4± would then further provide the matrix 
forms of Ai and A2, for example. 

A compact representation for the operators A is available within the degenerate 
subspace of dimension AN and corresponding to E"^ = N, which is the orthogonal sum 
of the two energy eigenspaces with the positive and negative energies of equal magnitude. 
It is straightforward to write down the operators H and J as block-diagonal matrices, 

H , , 

= ^3 ® , (27) 

V A* 

J, =h®U. (28) 

The identity matrix I2N and the projected (reducible) angular momentum matrices L 
act within the degenerate subspace of ii^ = \/iV, and oiE = —\fN. The chiral symmetry 
operator F that has the effect of reversing the sign of energy can be similarly expressed 
as 

T = ai®hN. (29) 

The vector operator A projected onto the E^ = N subspace may be seen to possess the 
form, 

2 

Ai = ^(73 ^ L, + 2VN(r2 ^ Ki , (30) 

in which through Eq. ( [T3|) the ffist term generates the above J. The second term 
contains yet another vector K, which due to the multiplication by a2 cancels out in 
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Eq. fll3p . but will be important for the analogy with the Runge-Lenz-Pauli vector to be 
drawn shortly. Via Eq. f|T9|) we then find 

[Ki, Lj] = ieijkKk , (31) 

which resembles the commutation relation between the components of the Runge-Lenz- 
Pauli vector and the orbital angular momentum in the hydrogen atom. [2] Moreover, the 
coincidence between and Af implies the following identity, 

-^2^" + ^r= ' (32) 

to be true for all i = 1,2,3. Consequently, the sum of squares for both vectors turns 
into L'^/N'^ + = 3/4. In the hydrogen's bound spectrum, on the other hand, the 
corresponding identity reads K"^ + L'^/n'^ = 1 — l/n^,[2] where the principle quantum 
number is denoted by n. The two identities indeed match, but only when the respective 
degenerate subspaces contain tvuo angular momentum multiplets. In hydrogen this is 
only true for the energy eigenspace with n = 2, which is 2s © 2p. In addition, one can 
also rewrite the anticommutators in Eq. (|T2|) in terms of the vectors L and K as. 



{A, A,} = h® [^^{U, L,} + 4N{Ki, if,}) , (33) 

from which the pair of anticommutators adds up to zero for different indices. For 
example, {Ki,K2} = (— i){[if2, -Z^s], -^^2} = {i/N'^){[L2, L3], L2} is obtained through 
utilizing Eq. (!3T|) and ( |32|) . and therefore the above sum of anticommutators vanishes. 



4. Effects of chemical potential and Zeeman field 

In this section we study the effects of chemical potential and of the Zeeman coupling 
to the external magnetic field perpendicular to the surface of the topological insulator. 
At the surface of topological insulator, the spin-rotational invariance is broken by the 
kinetic energy term in Eq. which suggests that the Hermitian matrix iaia2 should 
correspond to the out-of-plane component of the spin. Since the two mass-terms are 
rotated into each other by the matrix i/3i/32, we recognize this operator as representing 
the particle number, i. e. the generator of the U{1) symmetry associated with the 
superconducting phase. Therefore, the out-of-plane Zeeman field h and chemical 
potential fi shall enter the Hamiltonian through the couplings to iaia2 and i(3i/32, 
respectively. It will be useful, however, to represent these two terms by their sum 
\E' and their difference 11 instead: 

^ = ^{PiP2 + ^1^2) = i{a\a2 - 4ai) , (34) 

and 

n = ^{(3iP2 - aia2) = i{a\al + 0102) . (35) 
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We notice that the operator \1/ coincides with the fermion part of angular momentum 
operator J2 in Eq. f|T^ . so J2] = H"^] = 0. The operator 11, on the other hand, is 
a fermion "pairing" term, which is also SU{2) invariant: 

[n, j] = . (36) 

The chemical potential and the Zeeman coupling then can be written together as 
eU + with e = fi — h and 6 = fi + h. Before proceeding further, we exploit the SU (2) 
invariance of the Hamiltonian in Eq. (|2]) and rotate for convenience the operator \E' into 

\1/ = a|ai — a2'^2 , (37) 

which is the fermion piece of J3. 

Since \1/ and 11 act on the states by changing only the fermion part of wave functions, 
they can be viewed as two linear and orthogonal operators within the fermion subspace 
{|0),al|0),4|0),al4|0)}, in which the common, but for the present purposes irrelevant 
boson part is not shown explicitly. 11 projects onto the subspace {|0), ajaalO)} consisting 
of even number of fermions, since obviously naL2)|0) = 0- Moreover, 



n|0) = ia\al\0), (38) 
nal4|0) = -i\0) , (39) 

so n acts like the ct2 within this subspace. in contrast, is a projection operator onto 
the subspace {allO), agjO)} consisting of only one fermion: 

Wi|0) = 4|0), (40) 

^4|0) = -4|0), (41) 

\l/ therefore acts analogously to in this subspace. 

We may also introduce the projection operators from the chiral symmetry operator 
r, namely, 

p± = (i±r)/2 

because F anticommutes with both fermion operators ai^2- This can be seen, for example, 
from the fact that V-^-a\\0) = and V-a\\0) = a\\0) since the zero-mode has r|0) = |0). 
It follows that = UV- = 0, which can also be seen from the direct computation 
that yields 

= V+, ^2 = (42) 

The full Hamiltonian including the perturbations can therefore be written as, 

H' = H + V+{eIl)V+ + P_(5^)P_ . (43) 

It is easier to consider the effects of 11 and \I' separately. Since it commutes with J3 
and H^, the state \E,j,m) is by \1/ coupled only to the state with the same m and the 
same E"^ . As far as the coupling matrix for ^! in Eq. (1371) is concerned, the only relevant 
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'> '> > > > 



j^,m> 




Figure 1. The splitting of the degenerate multiplets |£',j>,m) and |i?,j'<,m) due to 
the effects of ^ in Eq. (gS]), for iV > 1 and fi ^ h = S/2. The thick sohd hnes stand for 
the two states with the same quantum number m but distinct j, which are coincident 
in energy for (5 = 0. The thin solid lines denote the unique states of maximal value of 
|?7i| within this degenerate subspace. Dashed line are the shifted energies when S 0. 
The middle dotted line stands for the unperturbed zero energy state. 



states are |1(3)) = \ ±E,j^,m) and |2(4)) = \ ±E,j^,m) when \m\ < j>. On the other 
hand, in the case of largest |m| = j>, the relevant states only include |1) and |3). It will 
be shown shortly that H + 5\1' is in this basis reduced to the universal matrix, 

Ecxs + 5{l2 - (Ti) ® (T3 , (44) 

when \m\ < j>. The two-dimensional matrix, 

Ea^ ± 5{h - ^1) , (45) 

on the other hand, is responsible for the cases of m = ±j>. In either case, the universal 
matrix is independent of m, and the eigenvalues are -^\/ E"^ + 5^ ± 5. To the first order 
of 5, the splitting of the spectrum's degeneracy is shown in FigJH 

Let us explain the appearance of the above universal matrix in more detail. First 
consider \m\ < j>. Recall that ^ acts only on the one-fermion part of wave functions. 
Hence, only the projected wave functions, 

P_| ±E,j,m) = ± [ci(2)/™(6l,4)ai±C2(i)^7™(&L4)4] |0) , (46) 

are needed when considering the effective couplings. The above and are some 
general products of boson operators 6^'s that will not affect the following analysis. 
Within the bracket, the order of coefficients Ci and C2 as well as ± sign is to denote 
the two orthogonal states corresponding to j = j> and j = j<, respectively, whereas 
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the ± sign preceding the bracket is to denote the sign of E. Now it is straightforward 
to show that \1/ is given by, 

S^^f = 6{l2 - , (47) 

in which the minus sign preceding cxi appears because opposite signs are associated with 
opposite-energy states in Eq. P6|) . The 2-by-2 matrix S is associated with couphng 
between |1) and |2). Exphcitly, 

S={cl- cDas + 2ciC2ai , (48) 

which yields a pair of eigenvalues ±{cl + cl). We note that + = 1 is dictated by the 
orthogonality between wave functions in Eq. f H6l) . Moreover, we can rotate the space 
spanned by |1) and |2) so that S turns into a diagonal cxs, which finally leads to the 
matrix form in Eq. ( l44l) . Similar method can be applied to \m\ = j> to obtain Eq. (I45l) . 

Next we consider the coupling matrix for 11 in Eq. fl35l) and the ensuing splitting 
of the spectrum. Since [H, j] =0 and [11, H] ^ 0, it suffices to consider another 
subspace consisting of |1(2)) = | ± Ei,j,m) and |3(4)) = | ± E2,j,m) for a given 
j with the constraint Ef = E^ + 2. Note that 11 will only act on the even-fermion 
part of wave functions P+|-E'i(2), j<(>), m). Explicitly, the one including both fermions 
aJa2Mm(&L ^2)10) corresponds to {E,j) = (i?i,j<), while the one without any fermion 
Mm(&i,^2)|0) is responsible when {E,j) = (_E'2,j>). As can be seen in the Appendix, 
Um represents a general product of boson operators 6^'s, which will neither affect the 
coupling matrix of 11. Then, the four-dimensional matrix representing H + eU assumes 
the following universal form, 

Esh 8)^3 + Ed(T3 ® + e (J2 ® (/2 + di) , (49) 

where the energy is defined through -Es(d) = (-^'i =t E2)/2. This matrix then yield two 
new pairs of eigenvalues, \E[\ = |ii^i|^/T+~|7p and |i?2| = \E2\^1 — I^P, respectively. 
As shown in Fig. |2l the degeneracy between a pair of distinct j is removed, while the 
accidental SU{2) degeneracy remains unresolved. 

The effects of the perturbations under consideration on the zero-mode can also be 
simply analyzed by using the projection operators. Obviously, due to the projection into 
odd-fermion subspace, the term proportional to in Eq. f H3|) by itself has no effect 
on the zero-mode. The 11 term in Eq. fH3l) can only couple the zero-mode to another 
pair of rotationally invariant states, i.e. the pair of opposite-energy states of = 2 
and j = 0. When the Il-term alone is present, the corresponding three-dimensional 
matrix that connects the j = states is then found to always have a zero eigenvalue. 
Within the linear approximation, we therefore conclude, the zero-mode is still present 
in the spectrum in the presence of the chemical potential and the Zeeman field, at the 
special point = \h\. This agrees with the more general analysis, [23] which shows 
that in the presence of chiral-symmetry breaking terms such as h and /i, the zero-mode 
is implied under very broad conditions by the existence of an antilinear operator that 
anticommutes with the Dirac vortex Hamiltonian. 
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|j =j.m> 



|j =j,m> 



E2(1-|er) 



E,(1+|er) 



Figure 2. The splitting between the states with the same quantum namber j but with 
different E, due to the Il-term in Eq. (|43|. This corresponds to ^ = —h = e/2. The 
sohd and dashed hues represent the multiplets before and after, respectively, the extra 
term e ^ is turned on. The dotted line in the middle stands for the unperturbed 
zero-energy state. The degeneracy due to the SU (2) rotational symmetry remains 
intact in this case. 



In presence of both terms \1/ and 11, the calculation becomes rather complicated, 
and we will not pursue it further here. As long as the perturbations have weak spatial 
variation, however, our treatment in terms of fermion operators only is still applicable. 
Presence of both perturbation in Eq. ( H3l) introduces couplings which are not closed 
within any finite subset of bound states. For each value of the quantum number m, 
there will be a corresponding matrix of infinite dimension, representing the coupling 
among all the states with the same m. In this case no degeneracy is expected to survive 
in the spectrum. 

5. Discussions and Conclusions 

The accidental symmetry and the emergent supersymmetry due to J and A both 
depend crucially on the assumed linearity as well as on isotropy of the vortex 
amplitude. Including the higher order terms in |A(r)| would introduce combinations 
of higher powers of boson operators, which would be expected to remove the accidental 
degeneracy. The anisotropy in the linear terms, on the other hand, would result in an 
asymmetry to the square of Hamiltonian H"^ = J2i=i,2 ^iio-ltti + bjbi), [19] and mi 7^ m2. 
Following the arguments in Sec. [21 it is evident that finding out the square roots of 
other than H itself is then impossible. Basically, it is due to the fact that the matrix 
(mi + 777,2 ) /2/2-I- (mi— 7712) cr3/2 has only one independent square root when the mi 7^ m2. 
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Terms in the Dirac Hamiltonian with the hnear dependence on the coordinates may 
also occur in graphene, for example, in presence of the magnetic field perpendicular to 
the plane. In this case a rather different supersymmetry was found. [2^ Yet another 
similar but maybe a better known case is the Dirac oscillator [25] with the Hamiltonian 

Has = ai{pi - iTxi). (50) 

In this case the crucial difference from our problem lies in that the four Dirac matrices 
{ai,a2,iaiT,ia2T} that appear in Hos do not close a Clifford algebra, i.e. not all the 
pairs of matrices anticommute. The square of the Hamiltonian, H^^, contains a spin- 
orbit coupling (1 — L35'3)r in addition to the common part of oscillator Hamiltonian as 
in our H'^. Here, the angular momentum L3 = xiP2 — X2P1 and spin S3 = iaia2 are 
defined. Compared to the present problem, the extra term in H^^ can be any integer 
and half-integer, whereas the extra term Yl,i=i,2 cf\cbi problem can have the values 

{0, 1, 2} only. These subtle differences ultimately lead to a rather different underlying 
(super) symmetry and the ensuing degeneracies of the spectrum. 

By solving the differential equation for the zero-mode's wave-function of our 
linearized Hamiltonian one finds its radial dependence to be exp (— AooT^/2c^). Thus, 
the linearization of the vortex profile is a good approximation when the coherence length 
is long, c -C Aoo^- Also, in the case of strictly linear vortex amphtude, strong external 
Zeeman field is not able to remove the zero-mode from the spectrum. This stands in 
contrast to the previously studied case of an order parameter which saturates to a finite 
value far from the vortex. [231 [26] The spinor part of zero-mode can be deduced from 
the matrix equation ai|0) = 02 10) = 0. Using the previous notations in Ref. [26], the 
fermion operators have the explicit matrix form, 

a. = { ? I a. = f \ . (51) 



2 



which yields the column vector |0) oc (1, 0, 0, 1)"^ for the zero-mode. 

Lastly, the treatment of chemical potential and Zeeman field with two projections 
\1/ and H is also relevant to a recent study on two-velocity Weyl fermions realized in 
two-dimensional optical lattice. [27] The effective low-energy Hamiltonian for the spinless 
fermions can be written as 

Hw = Kai + kya2 + rjiikJ^Pi + ikyT(32) , (52) 

which yields the eigenvalues ±(1 ± r])k, and r/ is a number indicating the difference 
in "speed of light". Actually the above Weyl Hamiltonian can be expressed as a 
more compact form as one notices that the sets of matrices U = {ai,a2,ic(i<y2} and 
V = {iTf3i,iTf32,i/3if32} form two commutating su{2) Lie algebras. Thus a generic 
Hamiltonian of 

Hg = u-U + V-V , (53) 

with components u = {k^, ky, 0) and v = rjik^, ky, 0), can be co-rotated into l-ulf/sH- l-ul V3. 
Then it is easy to see that Hw = (1 + i])k'$ + (1 — ri)kll, which explains that 
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the faster/slower Weyl fermions in fact correspond to the odd-fermion/even-fermion 
subspaces in our definition. 

In conclusion, for the two-dimensional linearized vortex on the superconducting 
surface of topological insulator, we defined the corresponding supersymmetric Runge- 
Lenz-Pauli operator, and discussed the connection to its counterpart in the hydrogen 
atom. We also studied the effects of the physical perturbations such as finite chemical 
potential and Zeeman coupling of the spin to the out-of plane magnetic field on the 
spectrum, and showed how these can be economically accounted for in terms of the 
even/odd- fermion projection operators which we introduced. In the most general case, 
the accidental degeneracy of the spectrum is removed by these perturbations, but the 
energy of the zero-mode is unaltered. Unexpectedly, however, we found that some of 
the accidental degeneracy remains present when the chemical potential and the Zeeman 
field have identical magnitudes. 
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6. Appendix 



Within the subspacc of states of E'^ = N, the wave functions with the maximum 
eigenvalue of J3 belong to the j>-multiplet. Explicitly, the pair of opposite-energy 
states can be written as. 



|0). 



(54) 



Within the same multiplet, the states with one lower in m is obtained by applying J_ 
to the above states. The even-fermion part of = ±y/N,jy,m — j> — 1) is. 



1 bM) 



= |0), 



1) 

while the associated odd-fermion part is, 

^W^alblibl)^-' 



(55) 



V2 



N 



+ 



1 alibi) 



N-2)\ \/N J{N-1)\ 



|0), 



(56) 



where the above ± denotes the corresponding sign of energy. In the other multiplet of 
_7<, the state \E = ±\/N,j^,m = j> — 1) has the same m as the above one, and they 
are orthogonal to each other. Therefore, its even-fermion part is given by. 



|0), 



(57) 
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while the associated odd-fermion part is, 

T -^^^=^—-\l^^ . |0), (5^ 



in which a reverse sign is used to denote the corresponding positive/negative energies 
in this multiplet. This choice is necessary for being consistent with the choice of sign 
for the unknown d in Eq. (1221) . One can find out the rest of wave functions within this 
subspace by successively applying the ladder operator J_ to the above states in both 
multiplets. 
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